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ABSTRACT 

We report on the detection of Main Sequence stars belonging to the recently identified 
Canis Major galaxy in a field located at ~ 4.2° from the center of the stellar system. 
With Main Sequence fitting we obtain a distance modulus (m — M) = 14.5 ± 0.3 
to the dwarf, corresponding to a distance of Dq ~ 8.0 ± 1.2 kpc, in full agreement 
with previous estimates based on the photometric parallax of M-giants. From the 
comparison with theoretical isochrones we constrain the age of the main population 
of the Canis Major system in the range ^4 — 10 Gyr. A blue plume of likely younger 
stars (age < 1 — 2 Gyr) is also identified. The available Colour Magnitude Diagrams 
of open clusters that may be projected onto the main body of Canis Major are also 
briefly analyzed. The position, distance and stellar population of the old open clusters 
AM-2 and Tombaugh 2 strongly suggest that they are physically associated with the 
Canis Major galaxy. Using our own photometry and data from 2MASS and the GSC2.2 
surveys we demonstrate that the claim by Momany et al. that the CMa overdensity 
is entirely due to the Galactic warp is not supported by the existing observations, 
once all the available pieces of information are taken into account. It is shown that 
the CMa overdensity clearly emerges at a heliocentric distance of ~ 8 kpc above any 
overdensity possibly produced by the Galactic warp. 

Key words: Galaxy: structure - galaxies: dwarf - open clusters: general - open 
clusters: individual: NGC 2477, Tombaugh 1, Tombaugh 2, AM-2, Berkeley 33 



1 INTRODUCTION 



The process of hierar chical merging JWhite fc ReesHl978t 
IWhite fc Frenk I Il99 ill is generally accepted as the driv- 
ing mechanism of the formation of giant galaxies. The 
study of the local (e.g. Galactic or in the Local Group) 
relics of such processes provides an unprecedentedly de- 
tailed insight of the current cosmological model as well as 
a formidable testbed for theories of galaxy formation. The 
case of the Sagittarius dwarf galaxy lllbata et alll994Lll997t 
Sgr dSph) is emblematic, since we are presently observing 
it during its disruption into a giant stream contributing 
to t he continued assembly of the Galactic halo with stars 
(see llbata et al.ll2001al : llbata et al. Il2002l llvezic et alJl200ti 



* Based on archive data from observations made with the ESO- 
MPI WFI telescope, at ESO, La Silla, Chile under programme 
170.A-0789(A) 
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iNewberg et al 1 l2002t iMaiewski etHI 12001 and references 
therein ) and globular clusters (jBeUaazini^^errar^^^bat 
iBellazzini et al.l l2003bl: Bellazzini. Ibata fc Ferrar 



2003a 



2003c ). A similar case has been reported in M3l1ibata^t^rt" 



2001b; iMcConnachie et al.ll2003T) and systematic searches 
in galaxies outsid e the Local Group are beginning 
jPohlen et al.ll2003lb 

An even more interesting case is provided by the 



recently discovered ring structure (the Monoceros Ring 



hereafter the Ring, for brevity jNewberj^j^ilJ 
[ Yannv et al. Il2003t llbata et ahllioOSl: iMaiewski et alJ 
iRocha-Pinto et alJ 



2002 



2003 



IRocha-Pinto et alJ I2OO3I) near the Galactic plane. The 
more r ecent observational results dCran e 2003 ; Martin et all 
l2004al) suggests that the Ring is the stream remnant of an 
in-plane accretion that may be contributing t o the b uild- 
up of the thick disk (see also lHelmi et al.ll2003l : IXbadi et alJ 
12003ft . 

In particular iMartin et alJ <l2004al . hereafter Pap-I), 
while tracing the whole extent of the Ring using M-giants 
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from the 2MASS database 1 discovered a large overdensity of 
"Ring- like" M-giants in the Canis Major constellation. The 
elliptical shape, the overall structure of this overdensity and 
its co incidence with a noticeable groupin g of globular clus- 
ters l|Bellazzini. Ibata fc FerrarJ |2003(J) led iMartin et alJ 
(2004a) to the conclusion that it is the relic of the dwarf 
galaxy whose disruption generated the Ring. The discov- 
ered relic (hereafter Canis Major galaxy, CMa) is approxi- 
mately centered at galactic coordinates (l,b)~ (240°, —7°), 
it has a FWHM extension in the sky of ~ 12° in the lat- 
itude direction and it is located at ~ 7 kpc from the Sun 
(i.e., (m — M)o ~ 14.25), as obtain ed from the photometri c 
parallax to M-giants introduced bv iMaiewski et all J2003I) . 
The mass, luminosity and characteristic dimensions of Canis 
Major appear quite similar to those of the Sgr dSph. The 
similarity with Sgr extends to other relevant characteristics: 
the two objects seem to host a similar grouping of globular 
clusters, and the upper RGB of both galaxies is dominated 
by stars of similar colour (M-giants), which are present in 
similar numbers in both galaxies. The analogy between the 
two systems may provide a useful guideline in the study of 
the newly discove r ed gala xy. 

IMartin et all J2004al) obtained a near infrared (NIR) 
Colour Magnitude Diagram (CMD) of the discovered struc- 
ture (e.g., their Fig. 8) showing a Red Giant Branch (RGB) 
extending to (J — Ks)o ~ 1.3 and a pronounced Red Clump 
(RC) at K so > 12.5 and (J - K s )o ^ 0.6. A Blue Plume 
(BP) of stars (of more uncertain nature) is also detected 
around (J - K s )o ^ 0.3. The Main Sequence (MS) of the 
CMa population is beyond the reach of 2MASS photome- 
try. The degeneracy among distance, age, metallicity and 
reddening (which is particularly high and highly variable 
at such low latitudes) prevents any firm conclusion about 
the nature of the stellar population in the CMa structure. 
Furthermore the distance scale based on p hotometric par- 
allax of M-giants is intrinsica lly uncertain feellazzini et alJ 
l2003bl : IMaiewski et alJl2003ft and it is likely affected by a 
system atic effect leading to the underestimate of true dis- 
tances iMartin et al.ll2004al: fNewberg et al. 1 12003ft . Hence, 
the detection of fainter stars belonging to the newly dis- 
covered galaxy, the access to alternative distance indicators 
than M-giants and the study of its stellar population in other 
wavelength ranges (than the NIR) are highly desirable to 
confirm the discovery presented in Pap-I and to obtain a 
clearer characterization of the system. 

Here we present the first detection of the Main Sequence 
(MS) of the CMa system in a wide field located ~ 4.2° apart 
from its center, from deep B and V photometry of very re- 
cent ESO-WFI images taken from the ESO archive. We also 
take advantage of the fact that the low latitude region of 
sky covered by the CMa galaxy hosts a number of galactic 
open clusters (OC), most of them in the foreground. Thus 
we searched for existing wide field CCD photometry of these 
open clusters to look for evidence of a background popula- 
tion that may be related to the CMa galaxy. Unfortunately, 
modern wide field photometries of open clusters are quite 
rare since, typically, a field of few arcmin is sufficient to 



1 See Cutri et al. (2003), Explanatory Supplc- 
ment to the 2MASS All Sky Data Re lease , 
http: / / www.ipac.caltech.edu/2mass /releases / allsky / doc / explsup. 



characterize the population of these clusters and the back- 
ground population is an undesired contaminant for these 
kind of studies. Nevertheless, some interesting clues of the 
presence of CMa stars in these fields have been obtained, 
providing useful indications on the spatial extent of the sys- 
tem and guidance for future observations. 

The list of the fields considered in the current study 
is reported in Table 1, where we report the position of the 
fields in Galactic coordinates, the field of view covered by the 
available observations, the angular distance from the center 
of CMa, an estimate of the reddening and the source of the 
data. The position of the various fields with respect to the 
CMa system is sketched in Fig. 1. Most of the photometry 
(catalogues in electronic form) used in th is paper have bee n 
retrieved from the WEBDA 2 database jMermirliodlll995ft . 
The adopted reddening value s are t aken from the maps by 
ISchleeel. Finkbeiner fc Davis! dl99ct SFD98), if not other- 
wi se stated. All over the paper we adopt the reddening laws 
bvlDean 3 _^^r£eri_fc_Cc^in^j [l978h for the (Cousins') I band 
and by |Rleke^^ebofsk^ ' lll985ft for all the other passbands. 

We anticipate that the data presented in this paper 
are not sufficient to definitely break the above quoted dis- 
tance/ 'age /metallicity /reddening degeneracy and that reli- 
able spectroscopic abundances and dedicated optical pho- 
tometry are still needed to obtain a fully consistent picture 
of the main characteristics of the galaxy. On the other hand, 
the results presented here should be regarded as a significant 
progress in the understanding of the CMa system as well as 
a confirmation of its existence. 

Criticizing the results presented in lMartin et al] <l2004aft 
and in a previous prelimin ary version of the pr esent pa- 
per circulated as a preprint. iMomanv et al] ]2004l . hereafter 
M04) claimed that the CMa overdensity is completely due to 
the North-South asymmetry produced by the Ga lactic warp 
(see lYusifovlliool iKuiiken fc Garcia-Ruizll200ll and refer- 
ences therein). We discuss this point in Sect. 4, demonstrat- 
ing that the M04 hypothesis is clearly confuted by several 
independent observational evidences. 



2 THE MS OF THE CANIS MAJOR GALAXY 

AT (L; B) = (244.2°; -8.2°). 

2.1 Observational material and data reduction 

We retrieved from the ESO archive a set of deep B and V 
images of a field located at (l;b)=(244.2°; -8.2°), at ~ 4.2° 
from the center of CMa. The images were taken with the 
WFI camera 3 mounted at the 2.2m telescope at ESO (La 
Silla, Chile) as part of the large programme Public Imag- 
ing Survey: WFI Follow-up of XMM-Newton Serendipitous 
Fields (170.A-0789(A)). We will refer to this field as the 
XMM Field (F-XMM), following the name of the origi- 
nal programme. All the calibrating images (bias, flat-fields, 



2 |http : //obsww . uni ge . ch/webda/webda . h tml | 

3 The WFI camera is a mosaic of eight 2048 X 4096 
pixels chips, with pixel scale 0.238 arsec/px and 
with a total field of view of ~ 34' X 33'. See 
http://www.ls.eso.org/lasilla/Telescopcs/2p2T/E2p2M/WFI/ 
for further details. 
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Table 1. Explored Fields. 



Field Name 


1° 


b° 


f.o.v. 


dist a 


E(B-V) 6 


Ref. 


F-XMM 


244.2 


-8.2 


34' 


x 33' 


4.2° 


0.14 


This work 


NGC 2477 (B,V) 


253.6 


-5.8 


34' 


x 33' 


13.6° 


0.77 


Momanv et al. (2001) c 


Tombaugh 1 


232.3 


-7.3 


7.0' 


x 7.0' 


7.7° 


0.43 


Carraro & Patat f 19951 


Berkeley 33 


225.4 


-4.6 


10.5' 


x 10.5' 


14.8° 


0.75 


Mazur, Kaluzny & Krzcmiiisky (1993) 


NGC 2204 


226.0 


-16.1 


14.7' 


x 14.7' 


15.9° 


0.10 


Kassis et al. (1997) 


Melotte 66 


259.6 


-14.3 


14.7' 


x 14.7' 


20.2° 


0.22 


Kassis et al. (1997) 


NGC 2243 


239.5 


-18.0 


15.0' 


x 15.0' 


10.5° 


0.07 


Kaluznv. Krzcmihskv & Mazur (1996) 




Open Clusters^ 














Arp-Madore 2 


248.1 


-5.9 


5.5' 


x 5.5' 


8.3° 


0.73 


Ortolani, Bica & Barbuv (1995) 


Tombaugh 2 


232.9 


-6.9 


10.4' 


x 10.4' 


6.9° 


0.39 


Kubiak et al. (1992) 



a Angular distance from the centroid of CMa, in degrees. 

^ From the reddening map s bvlSchlegel, Finkbeiner Or, Davisl ll998h 

c The (V,I) photometry by Kassis et al. 11997) has been also considered (f.o.v=14.7' X 14.7'). 
^ Possibly associated with CMa, see Sect. 5. 
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Figure 1. Relative positions on the sky (Galactic coordinates) of 
the CMa galaxy and of the fields considered in the present anal- 
ysis. The center of CMa is represented as a cross, the two con- 
centric circles have radius r= 6° deg and r= 12°, corresponding 
to the HWHM and FWHM of the distribution of CMa M-giants 
as derived in Pap-I, respectively. The clusters that are possible 
members of the CMa system are represented as open stars. 

standard star fields) related to the same observational run 
were also retrieved. 

To obtain stellar photometry of this field we selected 
the best-seeing, lowest airmass (AM < 1.1), V and B im- 
ages of the se t taken during a photometric night in which ob- 
servations of iLandoItl l|l992l) standard fields, taken with the 
same instrument, are also available (the night of February 1, 
2003). The seeing and the exposure time are 0.88" FWHM 
and 440 s for the V image, and 1.07" FWHM and 360 s 
for the B image, respectively. The images were corrected 
for bias and flat-field using standard procedures within 



the IRAF environment mscred and were reduced using: 
the P SF-fitting software DoPHOT JSchecter. Mateo fc Sahal 
1992). The images were searched for sources adopting a 
3cr threshold above the background and modeling the spa- 
tial variation of the PSF with a quadratic polynomial (es- 
sential l y the same approach adopted in iBellazzini et alJ 
l2002bl : IBellazzini. Ferraro fc Ibatal l200ol and described in 
more detail in those papers). Only sources classified as 
bona-fide stars by DoPHOT were retained in the final cat- 
alogue of instrumental magnitudes, that was obtained by 
cross-correlating the individual B and V catalogues. All the 
sources having a photometric error (in B or in V) larger 
than 3 times the average error at their magnitude were also 
excluded from the catalogue. The instrumental photome- 
try has been corrected for the geometric distortion affecting 
the WFI camera (du e to the effect of scattered light, see 
iManfroid etai]|200ll) with a simple empirical formula de- 
rived by M. Irwin from the photometry of several Landolt's 
standard fields, that provides corrections very similar to the 
complex model provided by ESO 4 . The degree of crowding 
of the XMM Field is quite low, hence it has been easy to 
select a large number (> 50) of bright and isolated stars to 
obtain robust estimates of aperture corrections. 

The instrumental magnitudes have been transformed to 
the Johnson- Cousins system with the calibrating equations 

V - v = -0.134(6- v) + 23.97 (r.m.s. = 0.017 mag) (1) 

B - 6 = 0.337(6 - v) + 24.77 (r.m.s. = 0.023 mag) (2) 

(where capital letters indicate calibrated magnitudes and 
lower case letters indicate instrumental magnitudes), ob- 
tained from observations of Lan dolt's fiel d SA9 8. The ex- 
tended list of standard stars by Stetson (2000) has been 
adopted. The atmospheric extinction coefficients estimated 
on April 13, 2003 have been taken from the WFI homepage. 



4 Seewww.ls.eso.org/lasilla/Telescopes/2p2T/E2p2M/WFI/zeropoints/ 



4 M. Bellazzini et al. 




Figure 2. CMD of the ~ 0.5° X 0.5° XMM field, located at ~ 4.2° from the center of the CMa system. The error bars at B - V ~ 2.2 
are the average photometric errors as a function of magnitude. The dotted, dashed and continuous lines indicate the locus of the CMD 
where the joint Completeness Factor (CF) is 80%, 50% and 5%, respectively. The joint CF as a function of V magnitude (computed at 
B — V = 1.0) is shown in the inset. 



2.2 Artificial star experiments 

To characterize the completeness of the F-XMM sample de- 
scribed above we performed a set of artificial star exper- 
iments. The artificial stars were extracted from a Lumi- 
nosity Function (LF) similar to the observed one, modeled 
with the same PSF model used for the photometry (includ- 
ing spatial variations), and added onto the original frames 
( ~ 2000 at each run and adopting the strategy described 
in lBellazzini et all ll2002allb|) . to avoid interferences between 
artificial stars) that were then reduced in the same way as 
the original scientific images. A total of ~ 30000 stars were 



simulated in the V and B images. The completeness fac- 
tor CF = N recovered / N a i mu i a ted as a function of magnitude 
was determined separately for each passband, the joint com- 
pleteness factor at any given colour is then obtained by the 
product CF = CF V x CF B . 



2.3 The Colour Magnitude Diagram 

The CMD of the XMM Field is shown in Fig. 2. The dotted, 
dashed and continuous lines indicate the locus of the CMD 
where the joint CF is 80%, 50% and 5%, respectively. In the 
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Figure 3. Compa rison with the observed CMD of F-XMM (upper left panel) and the synthetic CMD from the Galactic model by 
iRobin et alJ J2003I) . Upper right panel: raw synthetic CMD; lower left panel: synthetic CMD corrected for incompleteness; lower right 
panel: synthetic CMD corrected for incompleteness and normalized to have the same number of stars as the observed CMD in the 
selection box plotted in CMDs on the left-hand panels. 



following analysis we will ignore the sources fainter than the 
CF=5% limit. 

Even a first glance to Fig. 2 shows that the F- 
XMM CMD is dominated by an unexpected sequence run- 
ning from (V;B - V) ~ (18.0; 0.6) to (V; B - V) ~ 
(23.0; 1.3), superposed on the typical field population. The 
morphology of the feature is typical of the MS of an 
intermediate-old stellar population. The possible Turn Off 
point (TO) of this sequence, at V ~ 19 is slightly bluer 
than the blue edge of the vertical plume of TO field 
stars located at various distances along the line of sight 
jGilmore fc Wysd Il987t lUnavane. Wvse fc Gilmorel 1 19961 
iPrandoni et al.lll999l : iMorrison et alJl200fJl . that is visible 



for V < 18 around B — V ~ 0.45 in this diagram. Note that 
the same characteristic (i.e. TO stars bluer than the average 
field TO stars) was one of the key elements that allowed the 
discovery of the Ring structure, that was first identified as 
an excess of stars o f spectral type A-F llNewberg et al .120021 ; 
lYannv et al. II2OO3T) . The sequence (as well as the overall ap- 
pearance of_tlie_CMD)jsalso very similar to the Ring MS 
found bv llbata et alJ ll2005t) . As in the CMDs presented by 
these authors, the Sub Giant Branch (SGB) and the base of 
the Red Giant Branch (RGB) of the unexpected population 
are hidden by the contamination of Galactic field stars. The 
saturation level of the photometry, occurring at V ~ 15.0, 
prevents the detection of the Red Clump (RC) of Helium 
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burning stars associated with the population, that are likely 
brighter than this limit. 

Other noticeable features of the CMD are (1) a sparse 
blue plume at B - V < 0.4 and V < 18.5, and (2) the 
vertical wedge-shaped plume of the local M dwarfs around 
B - V ~ 1.6. 

2.4 Comparison with a Milky Way model 

Lacking a suitable Control Field to compare with F-XMM 
we are forced to rely on a synthetic model of the stellar 
populations of the Milky Way. We adopt, her e and in the 
following, the Galactic model by iRobin et alJ (120031 . here- 
after R03), from which it is possible to obtain synthetic 
CMDs for a given direction and field of view through a 
web interface. The model includes all the known Galactic 
components (thin and thick disc, halo and bulge) and takes 
into account the effects of the disc warp and flare. The last 
characteristic is particularly useful in the present context, 
since the (stellar) Galactic warp reaches its maximum pro- 
jected density at I ~ 270° in the southern Galactic hemi- 
sphere, n ot so far away from the patch of sky we are con- 
sider ing jDiorgovski fc Sosml Il989l : lLopez-Corredoira et alJ 
120021) . We have checked the warp model adopted by R03 
and we found it in excellent agreement with th e observations 
presen ted by lLopez-Corredoira et alJ (l2002h . IRobin et alJ 
(2003) tested their model against observations in several 
directions and suggest that it provides reliable predictions 
down to V ~ 22. In the R03 model the extinction is mod- 
eled by a diffuse thin disc. In the comparisons with observed 
data presented in the following we adjusted the reddening 
assigned to synthetic stars by the model to match the ob- 
served degree of extinction toward the considered direction. 

In Fig. 3 we compare the observed CMD of F-XMM 
with the predictions of the R03 models for a field of the same 
area and in the same direction. The comparison is limited to 
the range in which CF^ 50% in the observed sample, i.e. to 
V ^ 21.5. We included the effect of photometric "noise" by 
adding to each synthetic star a Smag extracted at random 
from a Gaussian distribution having a equal to the average 
photometric error of the observed stars at the considered 
magnitude. 

Even from a first glance to the upper panels of Fig. 3 
it is quite clear that the above-described sequence that 
dominates the F-XMM CMD has no counterpart in the 
Galactic model. In the latter CMD the vertical sequence of 
field MS stars reaches V — 19.0 and then abruptly bends 
to the red, changing into a sparser sloped band with a 
sharp faint/blue edge. This is a well known general fea- 
ture of the CMDs of Galactic fields, observed by various 
aut hors and predicted by a l l Galactic models (see, for exam- 
plc Castcllani et al . 2002; iBellazzini. Ferraro &i Buonannd 
| l999al:lPrandoni et alJl999l : lMonaco et al.l2002tllbata et alJ 
2003). On the other hand the observed "unknown" sequence 
curves gently, always remaining to the blue of the synthetic 
sequence down to V = 21.5 (and beyond, not shown here), 
and remains densely clustered all over its length. The ob- 
served CMD shows also a clear excess of Blue Plume (BP; 
at 0.0 < B - V < 0.5 and 16.0 < V < 18.5) stars with 
respect to the synthetic CMD. 

The direct comparison with the raw synthetic CMDs is 
not quantitatively correct since it neglects the incomplete- 
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Figure 4. Luminosity Functions of the observed sample (filled 
circles and continuous line) and of the synthetic sample (open 
circles and dotted line). In the upper panel the raw synthetic 
sample is considered, in the lower panel the synthetic sample has 
been corrected for incompleteness and normalized as described in 
Sect. 2.4 (see also Fig. 3). 



ness of the observed sample. In fact, the synthetic CMD is 
100% complete at all magnitudes. The lower left panel of 
Fig. 3 shows the synthetic CMD once the incompleteness is 
taken into account. We have associated to each star of the 
raw synthetic sample the joint CF factor of the observed 
sample at its magnitude. Then we extracted at random a 
number in the range [0.0,1.0] and we compared it with the 
computed CF. If the random number was smaller than or 
equal to the characteristic CF of the synthetic star the star 
was retained, otherwise it was excluded from the sample. 
The overall appearance is now more similar to the observed 
CMD, in terms of number of stars in regions of the CMD 
dominated by bona fide field stars. However, the synthetic 
CMD still shows an excess of stars of a factor ~ 1.4 within 
the box superposed on the left-hand panels of Fig. 3, enclos- 
ing the brightest part of the sequence of field MS-TO stars. 
The synthetic M dwarfs brighter than V ~ 20.0 are over- 
abundant by the same factor with respect to their observed 
counterparts. These comparisons suggest that the model is 
over-predicting the number of field stars in this direction. 
To take into account also this possible effect we picked up 
stars at random from the completeness-corrected synthetic 
CMD until we obtained a sample having the same number 
of stars populating the selection box as the observed one. 
The CMD of the "normalized" sample is shown in the lower 
right panel of Fig. 3. This plot provides a realistic view of 
the contamination of the observed CMD by field stars. It 
may be appreciated that the highest degree of contamina- 
tion is around the TO of the "unknown" sequence, around 
B - V ~ 0.6 and for V < 19.0. Any eventual Sub Giant 
Branch of the unexpected population is completely hidden 
by field dwarfs having B — V ^ 0.8. These facts constitute 
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Figure 5. Results of the MS-fitting proced ure with isochrones of different metallicity. Upper panel: [M/H] = —0.66; lower panel: 
[M/H] = -0.35. The isochrones are from lPietrinferni et all <2004l) . 
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the main limitation in obtaining strict constraints on the 
age of the newly-identified population. 

In Fig. 4 we compare the LFs of the F-XMM field and 
of the synthetic sample for all stars having B — V ^ 1.1, to 
exclude the Galactic M dwarfs from the comparison. In the 
upper panel of Fig. 4 the observed LF is compared to the raw 
model predictions while in the lower panel the comparison is 
with the LF of the completeness-corrected and normalized 
synthetic sample. Fig. 4 shows that there is a significant 
excess of stars in the observed CMD corresponding to the 
unknown MS population, even if incompleteness and den- 
sity normalization are not taken into account. The excess 
of observed stars arise between V ~ 18.0 and V ~ 19.0 
and persist down to the limiting magnitude of the observed 
sample (not shown here). 

From the above comparisons we conclude that in the 
XMM field there is indeed a conspicuous stellar population 
that is not comprised in the R03 Galactic Model. We identify 
this population with the newly discovered Canis Major stel- 
lar system, since the considered field is well within its main 
body, as described in Pap-I. Hence we identify the sequence 
dominating the F-XMM CMD as the MS of the Canis Ma- 
jor galaxy, revealed for the first time by the observations 
presented here. 

2.5 Distance, metallicity and age 

Wetry to constrain the distance to the CMa system in 
the XMM field by fitting the observed MS with theo- 
retical isochrones. Here and in the following we will use 
isochrones taken f rom t he very recent and complete set of 
IPietrinferni et al.l <2004 hereafter P04 5 ). 

The fact that the RGB of the galaxy is dominated by 
M-giants strongly sugge sts a mean me t allicity larger than 
[M/H] ~ -1 feellazzini et alJ l2003bl : iRocha-Pinto et alJ 
l2003l: iMartin et al.ll2004al) . In~particular ICrand fcOOSt de- 
rived [Fe/H] = —0.4 from low resolution spectra of Ring 
M-giants in the vicinity of the main body of CMa. Hence 
we will adopt isochrones of metallicity [M/H\ = —0.66 and 
[M/H] = —0.35 for the present analysis. The results of the 
MS-fitting procedure are shown in Fig. 5. The upper panel 
shows the comparison with isochrones with [M/H] = —0.66, 
the lower panel shows the comparison with [M/H] — —0.35 
isochrones. 

In both cases the best fit is obtained assuming E(B- 
V)=0.10, in reasonable agreement with the SFD98 red- 
dening maps (E(B-V)=0.14). Several authors have sug- 
gested that SFD98 reddening maps may ove restimate the 
actual reddening by a factor ~ 1.14 — 1.35 JStanekl W9& 
Bonifacio. Monai fc Beers! l2000l : iDutra et all 120031: ISurnJ 
2003). If this possible correction is taken into account the 
agreement between our best-fit reddening and the estimates 
of the SFD98 becomes excellent. 

The best-fit distance modulus is (m — M)o = 14.4 ± 
0.1 and (m - M) = 14.7 ± 0.1 when the [M/H] = -0.66 
and [M/H] = —0.35 isochrones are adopted, respectively. 
The value (m — M)o = 14.4 corresponds to a heliocentric 
distance of Dq = 7.6 ± 0.4 kpc, in full agreement with the 

5 The isochrones can be retrieved from 
http://www.te.astro.it/BASTI/index.php 



estimate of Pap-I, within the uncertainties (-D0 ~ 7.1 kpc). 
Similarly, (m — M)q = 14.7 ±0.1 corresponds to Dq = 
8.7 ± 0.4 kpc, ~ 20% larger than the estimate of Pap-I but 
still fully compatible with it, given the ~ 30% systematic 
uncertainty affecting the M-gi ants distance scale adopted in 
that work fsee lMartin et alJ2004aHNewberg et al. l2003t for 
discussion) . 

No stringent constraint on metallicity can be derived 
from the comparisons shown in Fig. 5. We just note that the 
lower distance modulus obtained with the [M/H] — —0.66 
isochrones put the theoretical RC locus just above the sat- 
uration limit of our photometry, in better agreement with 
observations compared to the [M/H] — —0.35 MS-fitting so- 
lution. It has also to be recalled that theoretical isochrones 
predict the location of the Zero Age Horizontal Branch, i.e. 
the faint edge of real distributions of Helium Burning stars. 

A significant lower limit to the age of the main popu- 
lation can be obtained by comparison with isochrones with 
the bright end of the observed MS. From Fig. 5 we obtain 
an age > 6 Gyr for [M/H] = —0.66 or an age > 4 Gyr for 
[M/H] = —0.35. On the other hand the corresponding upper 
limit is made greatly uncertain by the high degree of field 
contamination in this portion of the CMD. Ages as large as 
10 Gyr cannot be excluded with the present data. Indepen- 
dent of the quoted uncertainties it should be also considered 
that (a) the CMa system probably hosts stars of various ages 
and mctalli cities , as most of the local dwarf galaxies actually 
do iMateolll99d) . and (b) the galaxy has a sizable dimen- 
sion along the line of sight {<J( m ~M) ~ 0.5 mag, according 
to Pap-I), hence the observed CMD sequences may appear 
wider and less well-defined than in smaller or more distant 
systems just because of this differential-distance effect. Fi- 
nally, to fully populate the observed Blue Plume, stars as 
young as ~ 1 Gyr should also be present in CMa, suggest- 
ing a recent (minor) burst of star formation, similar to what 
observed in the Sgr dSp h jBjjla zzini. Ferraro fc Buonannol 
ll999bl:lLavden fc Saraiedinil200Cft ~ 

From the above comparisons it appears clear that more 
complete photometry (reaching the Tip of the RGB) and 
a suitable (observed) control field to perform a robust field 
decontamination are mandatory to obtain a deeper insight 
into the physical characteristics and the stellar content of 
the CMa system. Nevertheless, we have provided here the 
first clear detection of the MS of the galaxy, an independent 
estimate of its distance, in good agreement with the results 
of Pap-I, as well as the first broad constraint on the age of the 
main stellar populations of the system. In the following we 
will use the results obtained in this section as guidelines for 
the analysis of the CMD of populations possibly associated 
with the CMa/Ring system in the background of Galactic 
open clusters. 



3 CMA IN THE BACKGROUND OF OPEN 
CLUSTERS? 

We searched the WEBDA database for published CCD pho- 
tometry of open clusters lying in the broad region covered 
by the CMa galaxy, i.e. 220° < I < 260° and -20° < 6 < 0°, 
according to Pap-I. We found a handful of possibly useful 
datasets (with sufficiently deep CCD data over a field of at 
least a few arcmin across), that we briefly analyze in the 
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following. The considered photometries, being optimized to 
image nearby clusters, are too shallow to provide firm detec- 
tions (and/or non- detections) of the CMa population. Un- 
certainties in the reddening may also affect the analysis of 
these fields. Nevertheless some interesting clues of the pres- 
ence of CMa stars can still be found and may provide useful 
hints for future observations. 



3.1 The background of NGC 2477 

The open cluster NGC 2477 is located at (l,b)= (1=253.6°; 
b= — 5.8°). Its age is ~ 1 Gyr, its heliocentric distance 
is ~ 1.3 kpc {Ka ssis et al.lll997f). an d the metallicity is 
[Fe/H] ~ -0.5 jFriel fc Janes! Il993h . Unfortunately, the 
field of NGC 2477 is strongly affected by interstellar extinc- 
tion ((E(B - V)) = 0.77 acc ording SFD98), whose actual 
amount is very uncertai n fsee lHartwick. Hesser fc McClurel 
Il972l: iKassis et al]ll997l) and varies significantly across the 
observed field (with a standard deviation of 0.07 mag, 
SFD98). This constitutes the major source of uncertainty 
for the results obta i ned in the following. 

iMomanv et alJ i200ll) recently presented the B,V pho- 
tometry of a wide field (~ 34' x 33') centered on the cluster. 
The CMD obtained bv lMomanv et al.1 (1200 ill is shown in the 
upper left panel of Fig. 6. The cluster MS is a narrow diag- 
onal band crossing the CMD from B - V ~ 0.5, V ~ 13.0 
to B - V ~ 2.0, V ~ 22.0. The cluster RC is clearly visible 
at B — V ~ 1.3, V — 12.5. The most striking feature of the 
diagram is the large number of stars that are located along a 
wide sequence paralleling the cluster MS from B — V ~ 1.0, 
V ~ 18.0 down to B - V ~ 1.5, V ~ 22.0. Indeed the 
feature is reminiscent of the upper MS - Turn Off (TO) 
region of an intermediate-old stellar population, like that 
described in the previous section and attributed to the CMa 
system. Moreover, from the bright head of this putative se- 
quence departs a Blue Plume of stars (0.5 ^ (B — V) < 1.0, 
17.0 < V < 20.0) similar to what is observed in the XMM 
field described in the previous section. 

The upper right panel of Fig. 6 shows the corresponding 
synthetic CMD from the R03 model. The photometric un- 
certainties in the synthetic CMD have been implemented as 
described in Sect. 2.4, above. The diagonal line over-plotted 
on both the observed and synthetic CMDs is intended to 
separate the background MS population we are interested 
in from the bulk of the cluster stars. In the present case 
the synthetic CMD is quite similar to the observed one, if 
we neglect the sequences associated with the cluster, sug- 
gesting that the CMD is dominated by Galactic field pop- 
ulation. However we note that (a) the actual morphology 
of the main observed sequence is different from what pre- 
dicted by the model, (b) the BP population is much richer 
in the observed CMD, reaching also fainter magnitudes, (c) 
a sparse clump around V ~ 17.5 and B — V ~ 1.6 in the ob- 
served CMD seems to have no counterpart in the synthetic 
one (see below). 

In the lower panel of Fig. 6 the LF of the stars lying to 
the blue of the selection line are compared. A slight excess 
in the synthetic LF is apparent for V > 18.0, growing to- 
ward fainter magnitudes and reaching a factor ~ 1.4 at the 
limiting magnitude of the observed CMD. However, it has 
to be considered that while the synthetic sample is 100% 
complete at any magnitude the observed sample should be 
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Figure 6. Upper panels: observed CMD of the NGC 2477 field 
(left panel) and the corresponding synthetic CMD from the R03 
model (right panel). The selection line separates the background 
MS sequence form the bulk of the cluster stars. Lower panel: 
comparison between Luminosity Functions. The symbols are the 
same as in Fig. 4 except for the smaller open circles that represent 
the synthetic LF once corrected for the completeness function of 
F-XMM shifted by —1.0 mag to account for the different limiting 
magnitudes of the two samples. 



significantly incomplete for V > 20, i.e. within the last two 
magnitudes from the faint limit of the photometry. Under 
similar conditions (i.e. same instrument, similar seeing, sim- 
ilar exposure times) the completeness of the F-XMM sam- 
ple is always lower than 50% in the same range near the 
faint limit of the photometry. If we correct the synthetic LF 
for completeness (as done in Sect. 2.4, above) adopting the 
Completeness Function of F-XMM (shifted by —1.0 mag to 
partially account for the different limiting magnitudes, i.e. 
V ~ 22.0 for NGC 2477 and V ~ 23.5 for F-XMM), the 
synthetic LF (small open circles) falls well below the ob- 
served one for V > 19.5. Hence, the test concerning the MS 
populations is not conclusive. 

The case for the evolved population and for the Blue 
Plume is more interesting and it is briefly discussed below 
in comparison with the arguments provided by M04 to ex- 
plain the population in the background of NGC 2477 as 
completely due to the Galactic warp. M04 made the same 
comparison shown in Fig. 6 here but they selected from the 
synthetic catalogue only the stars with heliocentric distance 
between 7 and 9 kpc, claiming that all the sequences possibly 
attributable to CMa (including the BP and the RC features) 
"are fully reproduced by the Galactic warped model" (R03). 
The comparisons shown in Fig. 7 below demonstrates that 
this is not actually the case. 

The upper panels of Fig. 7 shows once again the com- 
parison between the observed CMD of NGC 2477 and 
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Figure 7. Upper left panel: observed CMD of NGC 2477. The 
heavy dots are the RGB and RC stars selected in the infrared 
CMD shown in the lower left panel. Upper right panel: synthetic 
CMD of an equivalent area corrected for extinction with the same 
assumptions adopted by M04. The heavy dots are the synthetic 
stars with heliocentric distance between 7 and 9 kpc (as in M04). 
Lower right panel: smoothed histograms of the stars in the color 
range of the Red Clump. 



its synthetic counterpart from R03 . We corrected the 
synthetic CMD with the same assumptions on reddening 
adopted by M04 (< E(B - V) >= 0.52 ± 0.04, after 
iBonifacio. Monai fc Beers) l2000h . In the observed CMD we 
plot as heavy dots the putative RC and RGB stars of CMa 
we have selected from the Near Infrared CMD of the same 
region that is shown in the lower left panel of Fig. 7. The 
RGB sequence is readily identified in this diagram and also 
the RC peak is more evident (see also Pap-I and references 
therein). The triangular box encloses the bluest part of the 
BP distribution (which extends to B — V ~ 1.0 in the ob- 
served CMD). In the synthetic CMD we plotted as heavy 
dots the stars fulfilling the same selection criteria adopted 
by M04 (7 < d e < 9 kpc). 

The overall morphology of the observed CMD is reason- 
ably reproduced by the distance-selected synthetic sample, 
though there are still significant differences: 

(i) The BP sequence is not correctly reproduced in po- 
sition and in number. The observed stars falling in the se- 
lection box are 582 while the synthetic stars are just 285, a 
full 10 a difference. Moreover the faint envelope of the syn- 
thetic BP sequence is one full magnitude brighter than the 
corresponding feature in the observed diagram. 



6 Actually a different realization with respect to that shown in 
Fig. 6. To check the effect of statistic noise on simulated CMD we 
analyzed four different realizations of the same synthetic sample 
in the present case. The differences among the various realizations 
are negligible under any relevant aspect. 



(ii) The selection of synthetic stars in heliocentric dis- 
tance indeed show evidence for a Red Clump population. 
However, as clearly shown by the LFs in the lower right 
panel of Fig. 7, the synthetic RC peak is more than 0.8 mag 
brighter than the observed one. Note that the positions of 
the peaks are unaffected by the inclusion of photometric er- 
rors in the synthetic sample. The synthetic RC appears more 
clumpy and well-defined in the plots produced by M04 (see 
their Fig. 3) because they neglect the effects of observational 
errors, but the average magnitude of the peak is the same as 
found here. Moreover, the overall comparison between the 
RC LFs suggests that the model is seriously overestimat- 
ing the total number of stars in the field (as for F-XMM, 
see Sect. 2.4), providing a further reason to reconsider the 
comparison of MS LFs shown in Fig. 6. 

(iii) M04 found that the average metallicity of the se- 
lected synthetic sample (heavy dots) is [Fe/H] = —0.45 ± 
0.25 and the average age is 5.0 ± 1.5 Gyr, similar to what is 
found for the CMa system. We verified these results with 
four different realizations of the synthetic sample of the 
NGC 2477 field and we fully confirm the metallicity result. 
However, it transpired that M04 misinterpreted the column 
of the output files of the R03 model labeled "age" as pro- 
viding the age in Gyr of any given star. In fact this is an in- 
teger code number that classifies the stars according to the 
Galactic component they belong to (thin disc, thick disc, 
bulge and halo) and to age ranges, within the thin disc pop- 
ulation (age/population class). Due to an incredibly lucky 
coincidence the average of the age/population class and of 
the actual age of the selected sample are the same (i.e. 5, 
a pure number, and 5 Gyr, respectively) but the standard 
deviations are significantly different (1.5 and 3 Gyr, respec- 
tively). A closer look at the age/population class distribu- 
tion shows that ~ 45% of the selected stars have age < 4 
Gyr, hence not similar to the typical CMa population, and 
~ 10% are thick disc stars, hence not possibly related to 
the warp. In conclusion, only less than 50% of the selected 
sample claimed to reproduce the CMa features can be real- 
istically related to it. 

(iv) Finally, we note that there is no reasonable reddening 
value that may simultaneously allow a reasonable match to 
all the main features of the synthetic and observed CMDs 
of Fig. 7 (i.e. MS, BP and RC). 

In conclusion the Galactic warped model adopted by 
M04 completely fails to reproduce the observations for what 
concerns the Blue Plume and the Red Clump in the back- 
ground of NGC 2477. The possibility that these stars are as- 
sociated with the CMa system or its debris (Ring) remains 
fully open. Deeper observations and extensive artificial stars 
tests are needed to verify the existence of their MS counter- 
parts at V > 18.0. 



3.2 The background of Tombaugh 1 and 
Berkeley 33 

The open cluster Tombaugh 1 ( Tom 1) is located at 
(l,b) = (232.2°, -7.3°). According to JCarraro fc PatatJll995l) 
the most likely metallicity is [M/H] ~ 0.0, the age is ~ 1 
Gyr and it is located at ~ 3 kpc from the Sun [(m — M)o ~ 
12.4]. In t he upper panels of Fig. 8 the CMD of Tom 1, as de- 
rived by llCarraro fc Patatill995|) . is compared to the corre- 
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Tom 1; (l,b) = (232.22,-7,32); Carraro & Patat (1995); FOV=7'x7' Be 33; (l,b) = (225.4,-4.6); Mazur et al. (1993); FOV= 10.5'x 10.5 1 
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Figure 8. Comparison of the observed and synthetic CMD (up- 
per panels) and LFs (lower panel) for the field of Tombaugh 1. 
The symbols are the same as in Fig. 6. 



Figure 9. Comparison of the observed and synthetic CMD (up- 
per panels) and LFs (lo wer panel) for the field of Berkeley 33 . 
The photometry is from iMazur. Kaluznv fc Krzeminskvl Jl993t) . 
The symbols are the same as in Fig. 6 and Fig. 8. 



sponding synthetic CMD from the R03 model. The observed 
CMD is dominated by the cluster MS going from V — I ~ 0.7 
and V ~ 14 to V-I ~ 2.0 and V ~ 21.5. Note that a broad 
and sparse parallel sequence seems to appear to the blue of 
the cluster MS (and of the selection line introduced in the 
previous section) at V ^ 18.5. The sequence seems to have 
no clear counterpart in the synthetic CMD and resembles 
the upper-MS of CMa as identified in F-XMM. The LFs 
shown in the lower panel of Fig. 8 confirms this view. The 
selected stars (i.e. the stars bluer than the selection line) 
are more numerous in the observed sample compared to the 
synthetic one all over the considered range of magnitudes 
(19.5 ^ V < 21.5). Also in this case the observed sample 
should be drastically pruned by incompleteness in the con- 
sidered range, so the real overdensity should be significantly 
larger than what is reported in Fig. 8. 

The same comparison is shown in Fig. 9 for the field 
of the open cluster Berkeley 33 (Be 33). Be 33 is located 
at (l,b) = (225.4°,-4.6°). It has an age of ~ 0.7 Gyr and 
[M/H] ~ —0.7; its distance from us is ~ 3.7 kpc (i.e ., 
\rn - M)o ~ 13.3, IMazur. Kaluznv fc Krzemiriskvl Il993l . 
The cluster MS crosses the diagrams from B — V — 0.6 
and V ~ 13 to B - V ~ 1.4 and V ~ 21. 

Also in this case a wedge-shaped sequence possibly cor- 
responding to the upper-MS of CMa is observed to the blue 
of the selection line and for V > 19.0. In this case the LFs 
of the observed and synthetic sample, shown in the lower 
panel of Fig. 9, are essentially indistinguishable. Again, the 
high degree of incompleteness that is expected to affect the 
observed sample for V ^ 20.0 suggests the presence of an 
excess of stars not predicted by the Galactic model. 

More complete and deeper photometry of wider fields 
surrounding the clusters are needed to confirm the tentative 
detections described above. 



3.3 Fields at b < -14° 

Finally, Fig. 10 shows the CMDs of three OCs hav- 
ing 220° < I < 260° and b < -14°, namely 
NGC 2204, Melotte 66 llKassis et al]ll9 97lL and NGC 2243 
llKaluznv. Krzeminskv fc M ! azu^ ~ Q99*6TL In spite of the 
larger field of view with respect to the observations of Tom 1 
and Be 33 described above, there is no sign of a defined 
upper-MS sequence to the blue of the selection lines in these 
CMDs. The assessment of the significance of a non- detection 
is a prohibitive task given the available observational mate- 
rial. The comparison with the CMDs of Tom 1 and Be 33 
(note that NGC 2204 is essentially at the same galactic lon- 
gitude as Be 33) simply suggests that the density of CMa 
stars is significantly decreased going from —7° ^ b ^ —4° 
to b < -14°. 

From the analysis presented in this section we simply 
conclude that the existing CMDs sampling the stellar pop- 
ulation in the background of clusters located in the same 
broad region of the sky that hosts the CMa system are com- 
patible with the results presented in Pap-I regarding the 
overall structure of the galaxy. As a guideline for future 
observations it is also suggested that the additional effort 
required to observe wider fields and to obtain deeper pho- 
tometry of clusters in this region may be worth doing and 
highly rewarding in terms of scientific returns, and it may 
also help to obtain a more thorough understanding of the 
population of the clusters themselves. 

4 THE CMA OVERDENSITY AND THE 
GALACTIC WARP 

Reconsidering the results of Martin et al. (2004a), M04 stud- 
ied the North/South asymmetry in starcounts of M-giants 
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Table 2. Selected Fields. 



Field Name 




1 


b 




r 


< E(B - V) > a 


Data Source 


CMa 


238° 


< I < 244° 


-11° < b < 


-6° 




0.27 ±0.09 


2MASS - PSC 6 


CMa-CF 


238° 


< I ^ 244° 


+6° < b < ■+ 


11° 




0.12 ±0.02 


2MASS - PSC 


WARP 


272° 


sC I «C 278° 


-11° < b 


-6° 




0.29 ±0.05 


2MASS - PSC 


WARP-CF 


272° 


«C I «C 278° 


+6° s: b sj 4 


11° 




0.23 ±0.04 


2MASS - PSC 


CMa (GSC2) 


i 


= 243° 


b = -8° 




r = 1.0° 


0.18 ±0.02 


GSC2.2 C 


CMa-CF (GSC2) 


i 


= 243° 


b = +8° 




T = 1.0° 


0.12 ±0.01 


GSC2.2 


WARP (GSC2) 


I 


= 274° 


6 = -8° 




r = 1.0° 


0.27 ±0.04 


GSC2.2 


WARP-CF (GSC2) 


I 


= 274° 


6 = +8° 




r = 1.0° 


0.21 ±0.05 


GSC2.2 



a Average reddening ± standard deviation from star-by-star interpolation on the SFD98 maps. 
^ Point Source Catalogue 

c Only entries classified as stars have been retained in all the samples extracted from the GSC2.2 catalogue. 



NG-C 2J304 
F0V=14.7'xl4.7 l 



Mel 86 

F0V=14.7'xl4.7' 



NGC 3843 
FDV=15'xl5' 




not the case with the only global model adopted by M04, i.e. 
the warped and flared Galactic model by R03 they use for 
comparison with the NGC 2477 field. Even accepting that 
this model reproduces the observations in the latter field 
(that is far from proven, see Sect. 3.1 above) it completely 
fails to reproduce the morphology of the CMD and the star 
counts in the XMM field, as exhaustively shown in Sect. 2.4, 
above. There, a very significant unexplained overdensity and 
stellar population remains, in excess of the warp model en- 
capsulated in the R03 model. Moreover in their analysis M04 
completely neglect the distance information provided by the 
photometric parallax of M-giants that is one of the key ele- 
ments of the discovery presented in Pap-I. In the following 
we shall demonstrate that the consideration of this addi- 
tional information easily singles out the CMa overdensity 
from the background of the Galactic warp. 

While extensively checked with various sets of observa- 
tions the R03 model cannot provide a completely safe com- 
parison and the ideal testbed to proof the reality of the CMa 
structure is the real Galaxy. In the following subsections we 
perform direct comparisons with suitable control fields using 
data from 2MASS and from the GSC2.2 8 surveys. 



Figure 10. CMDs of open clusters having 220° < I < 260° 
and b < -14°. The dat a for NGC 2204 and Melotte 66 are 
from iKassis et al.l Jl997fl. the data for NGC 2243 are from 
iKaluznv. Krzemihskv &; Mazurl Jl996t) . 



in the latitude stripe \b\ < 20° and in the latitude range 
235° ^ I ^ 245° and showed that in this spot of the sky the 
CMa overdensity is completely erased if (i) a model warp of 
amplitude w = 2.0° ±0.3° is adopted and (ii) th e correction 
introduced by iBonifacio. Monai fc Beers! J2Q00) is applied 
to the SFD98 reddening values. It is quite obvious that any 
overdensity in the proximity of the Galactic Plane can be 
locally reproduced by a particular model of the warp. What 
should be demonstrated is that a global model of the warp, 
broadly reproducing the observations in the whole Galaxy, is 
also able to reproduce the CMa overdensity. 7 In fact this is 

7 It should be noted that the warp model used by M04 is so crude 
that can be used only on small scale, by definition. To emulate the 
local effects of the Galactic warp they adopt b = —2° as the plane 



4.1 CMa is not the warp: 2MASS data 

All the tracers adopted in the literature to study the Galac- 
tic warp agree in placing the maximum amplitude of the 
southern Galactic warp around I = 270°, the North/South 
asymmetry remainin g near the max i mum in the broad range 
230° < I < 280° llYusifovl l2004t lLopez-Corredoira et alJ 
2002). The classical parameter adopted to quantify the 
asymmetry produced by the warp at a given (I; \b\) is the 
ratio between the counts (of any given density tracer) in 
the northern Galactic hemisphere and the counts in the 
corresponding region in the southern hemisphere (I;— \b\), 
R — N N or th / N south- At the longitudes of the maximum 
amplitude of the warp, R is ~ 0.7 in the southern warp 



of symm etry of the Galaxy, instead of b = 0° . In Martin et alJ 
12004b) we show that the CMa overdensity is detected even using 
the M04 warp model, once the distance information is properly 
taken into account. 

8 See http://www-gsss.stsci.edu/gsc/gsc2/GSC2home.htm 
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Figure 11. NIR CMDs of the 2MASS selected fields described in 
Tab. 2. The vertical lines encloses the color range populated by 
Red Clump stars at any distance. The adopted M giants selection 
box is plotted in the upper right panel. The arrow indicates that 
the box is shifted along the magnitude axis to count M giants as 
a function of Kq. 

jYusifovl Eool lLopez-Corredoira et alJ 120021) . Here, since 
we deal with southern structures we will use Rl = -g = 
Nsouth/NNorth- The typical value of Rl around the south- 
ern maximum is Rl ~ 1.4, which we will adopt as a reference 
value in the following. In any case values of Rl significantly 
larger than 1.0 indicate an asymmetry between the Galactic 
hemispheres, in particular an excess of star counts in the 
Southern hemisphere. Note that computing Rl from a ho- 
mogeneous dataset (as, for instance 2MASS) removes most 
of the problems associated with the incompleteness of the 
samples since the parameter is the ratio of starcounts from 
samples with very similar completeness as a function of mag- 
nitude. 

The key point of the following analysis is that if the 
CMa overdensity is due to the southern Galactic warp 
we shall observe similar overdensities and populations also 
around I = 270° where the warp asymmetry should be simi- 
lar if not larger. To do this comparison we carefully selected 
two windows of low reddening E(B — V) < 0.5 and of low 
average reddening {E(B — V)) < 0.3 and their respective 
control fields in the northern hemisphere with the same red- 
dening constraints. The data have been taken from 2MASS 
adopting the same selections as in lBellazzini et alJ (|2p03b). 
The chosen windows - of the same area and covering the 
same latitude ranges - are described in Table 2. For brevity 
the fields centered at I = 241° have been dubbed CMa 
(South) and CMa-CF (North), those centered at I = 275° 
are called WARP (South) and WARP-CF (North). The data 
have been corrected for extinction according to the SFD98 
maps. We note that all the resul ts presented in this section 
remai n valid also if the recipe bv lBonifacio. Monai fc Beera 
(2000) is adopted. The CMDs of the four considered fields 
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Figure 12. R/ distributions of RC stars in the CMa (filled circles) 
and WARP (open circles) fields in different ranges of Galactic 
latitudes. The dotted line marks the level of asymmetry typical 
of the Galactic warp, ij/1.4 ( Yusifov 2004; Lopez-Corredoira et al. 
120021) . 



are shown in Fig. 11. The vertical lines enclose the color 
range populated by RC stars at various distances from 
us. The parallelogram plotted in the uppe r right panel of 
Fig. 11 is the same selection box adopted by B ellazzini et alJ 
(2003b) to search for overdensities of M-giants as a function 
of distance. The box, shaped around the RGB of the Sgr 
dSph galaxy, is shifted along the magnitude axis counting 
the stars that fall into it as a function of the shift. The shifts 
are converted to di stance moduli assuming (m — M)p = 16.90 
for Sgr , as done bv lMaiewski et all (l2003Tl and lMartin et all 
ll2004al) . 

In Fig. 12 the distributions of Rl as a function of Ko 
magnitude and computed for stars in the color range of the 
Red Clump are plotted for different ranges of Galactic lat- 
itudes. In the range 6° \b\ ^ 7° a strong peak in Rl, 
reaching Rl — 2.1, is clearly detected in t he CMa direction 
at Kg ~ 12.7. According to the model bv lSalaris fc Girardil 
(2002), the absolute K magnitude of the Red Clump of a 
population having [M/H] ~ —0.4/ — 0.7 in the age range 
4-10 Gyr is < M K >=± -1.5 ± 0.2. Hence the peak ob- 
served in the CMa direction corresponds to an overdensity 
at the distance (m — M)o — 14.2 ± 0.4, in good agreement 
with the distance estimates for the CMa galaxy obtained in 
Pap-I and in Sect. 2.5, above. Note also that the Full Width 
at Half Maximum of the Rl peak is fully consistent with 
the FWH M of the CMa struc ture obtained in Pap-I (see 
Fig. 6, in lMartin et alJl2004ah . On the other hand in the 
WARP direction, while Rl > 1 essentially everywhere, as ex- 
pected, a much weaker peak is detected around Kq ~ 13.1. 
Let us assume that both the above described peaks are due 
to the Galactic warp. If we imagine the Galactic disc as a 
series of concentric circular rin gs with the o utermost rings 
tilted (because of the warp, see lBinnevlll992L and references 
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Figure 13. R/ differences from M giants as a function of distance 
modulus (see Bellazzini et al. (2003b) for details on the adopted 
technique). 

therein) the structures at / ~ 241° (CMa) and at I ~ 275° 
(WARP) should lie at the same Galactocentric distance, i.e. 
in the same ring. Seen from the Sun, this ring should appear 
more distant in the WARP direction, as observed, given the 
fainter RC peak detected in the WARP field. However, to 
obtain the observed difference in the magnitude of the peaks 
(AKo — 0.4 mag) the warped ring should lie at a Galacto- 
centric distance of ~ 20 kpc, implying a distance of the CMa 
structure from the Sun of ~ 15 kpc, in sharp contrast with 
the observations (see Pap-I and Sect. 2.5, above). On the 
other hand, assuming a realistic Sun-CMa distance of 8.0 
kpc, the difference between the two peaks should be ~ 0.9 
mag, much larger than observed. The only way to save the 
original hypothesis, that the overdensities in the CMa and 
WARP directions belong to the same warped ring of the 
Galactic disc, is to admit that such a ring is not circular, 
but significantly elliptical instead (axis ratio ~ 1.4). Even 
in this (very unlikely) case it would be very hard to explain 
the large difference in the amplitude of the Rl peaks, much 
higher toward I ~ 241° than toward I ~ 275°. 

Even more interesting, in this sense, is the behaviour 
of the Rl distributions at 8° < \b\ sC 9° (Fig 12, middle 
panel). Here there is no more peak in the WARP direction. 
For K > 12.0 the Rl distribution of the I ~ 275° field set- 
tles at the typical "warp" value (dotted line) while the CMa 
peak is as strong and well denned as at lower latitudes and 
in the same position. This comparison alone shows that (a) 
toward (/; 6) ~ (240°; -8.5°) the South/North asymmetry 
in RC stars is significantly larger than in the WARP field (a 
difference larger than la at Ko — 12.7) and (b) the excess 
of counts in the CMa direction is confined to a well denned 
range of magnitudes corresponding to the distance of the 
CMa galaxy as obtained in Pap-I and in Sect. 2.5, here. It 
should be concluded that while around I ~ 275° the distri- 
bution of RC stars is fully consistent with what is expected 



from the Galactic warp (see lLopez-Corredoira et al]l2002T) . 
in the I ~ 240° direction a clear and unexpected overdensity 
is observed at the distance of ~ 8 kpc, fully consistent with 
the results of Pap-I. It should be noted that the overdensity 
detected in the CMa direction is not only a South/North 
excess of starcounts in a limited range of distances but it is 
a clear excess with respect to what observed in the real warp, 
i.e. in the I ~ 275° direction. The result is confirmed also 
in the latitude range 9° ^ \b\ ^ 11°, where the CMa peak 
begins to weaken, as expected (see Fig. 6 of Pap-I) while the 
WARP distribution is essentially unchanged (in a greement 
with the results bv lLopez-Corredoira et al]l2002l see their 
Fig. 15). 

Finally, in the hypothesis that the South/North asym- 
metries observed in the CMa and WARP fields are both due 
to the Galactic warp, one can object that since the warp 
ring has a larger heliocentric distance in the I ~ 275° direc- 
tion than toward CMa, the same |6| range samples different 
ranges in Z (linear distance from the Galactic plane) in the 
two directions. It can be conceived that a comparison of the 
Rl distributions at the same Z is more appropriate than the 
comparisons shown in Fig. 12. Let us consider two possible 
scenarios: 

(i) If we assume that the warp lies on a circular ring at 
the same Galactocentric distance as the observed CMa over- 
density then the 9° ^ \b\ ^ 11° range toward CMa sam- 
ples the same range in Z than the 6° ^ \b\ ^7° range in 
the WARP direction. Comparing the CMa distribution in 
the lower panel of Fig. 12 with that of the WARP field in 
the upper panel of the same figure it can be noted that 
the amplitude of the peaks is similar (the CMa peak being 
still the strongest of the two) but the difference in magni- 
tude between the peaks is again 0.4 mag, which would force 
the conclusion that the warp is elliptical (already discussed 
above) . 

(ii) If we assume that the real difference in distance mod- 
ulus between the warp as seen toward I ~ 241° and to- 
ward I ~ 275° is 0.4 mag, as indicated by the position of 
the Rl peaks in the upper panel of Fig. 12, then the com- 
parison at the "same Z" should be done between the CMa 
distribution at 8° < |6| < 9° and the WARP distribution 
at 6° ^ \b\ ^ 7°. Since the CMa Rl distribution is essen- 
tially unchanged in the two latitude strips (upper and mid- 
dle panels of Fig. 12) this hypothesis must be rejected: the 
CMa peak is much higher and at a different Galactocentric 
distance with respect to the peak at I ~ 275° (possibly) at- 
tributable to the galactic warp. Note that this peak is seen 
only at 6° \b\ ^ 7°, while it is completely lacking at higher 
latitudes. 

Hence, the consideration of this possible effect is of no 
help to the M04 hypothesis: the observed Rl distribution of 
RC stars in the CMa field cannot be reconciled with the Rl 
distribution in the real warp (WARP field) . 

Essentially the same results are obtained from M giants, 
as shown in Fig. 13. In this case, since the tracers are much 
less numerous than RC stars and the signal-to-noise ratio 
of any structure is consequently lower, we had to adopt a 
coarser grid to compute Rl and the final result is best appre- 
ciated when shown as the difference between the Rl distribu- 
tions in the considered directions (Rt(CMa) - Rt(WARP)) 
vs. (m — M)o). It is clearly evident that the only significant 
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excess is the broad peak in the range 14.0 ^ (m— M)o 15.0 
that is detected in all of the considered latitude ranges. Note 
that Fig. 13 is directly showing excesses in star counts in the 
CMa direction with respect to the WARP direction, hence the 
CMa clump of stars clearly emerges over the South/North 
overdensity due to the warp. 

It must be concluded that the observed Rf distribu- 
tions of RC stars and M giants toward I ~ 241° and to- 
ward I ~ 275° are in clear contrast with the claim that the 
overdensity discovered in Pap-I and studied here is due to 
the Galactic warp (M04). On the other hand, the observa- 
tional evidence confirms that around (l;b) ~ (240°; —7.5°), 
at ~ 8.0 kpc from the Sun, there is a previously unknown 
stellar system that is clearly detected even above the under- 
lying asymmetry due to the Galactic warp: the Canis Major 
galaxy. 

4.2 CMa is not the warp: GSC2.2 data 

The optical (photographic) photometry provided by the 
GSC2.2 catalogue is a much less homogeneous dataset with 
respect to 2MASS. The effect of interstellar extinction is 
much stronger in the GSC2.2 passbands {Jmag,F mag ) and 
the corresponding reddening laws are not well known. On the 
other hand a direct test performed on the F-XMM showed 
that the limiting magnitude of GSC2.2 photometry is ~ 1.0 
mag fainter than that of 2MASS at B - V «S 1.0. Hence 
the inspection of the GSC2.2 data may provide additional 
insight in the comparison between the CMa overdensity and 
the Galactic warp. We have selected from the GSC2.2 two 
small circular regions of radius r = 1° toward the CMa 
and WARP directions and their respective northern Con- 
trol Fields. The fields have been selected to have small and 
homogeneous reddening (see Tab. 2, for details). Compari- 
son among the overall LF of the four fields indicates that 
the completeness of the samples is very similar down to 
Jmag ^ 17.5. Small color shifts (< 0.2 mag) have been 
applied to match the main features of the corresponding 
CMDs. We stress that the result discussed below is com- 
pletely independent of the assumed shifts. 

The CMDs of the selected fields are shown in Fig. 14. 
A detailed analysis of the diagrams suggests that there are 
differences in the stellar populations of the various fields, 
particularly at faint magnitudes, and the study of LFs of 
selected populations provides results fully compatible with 
those described in the previous section, above. However, 
given the possible problems with reddening corrections and 
the uncertain homogeneity of the samples 9 we drop this part 
of the analysis and concentrate on a single feature of the 
CMD, i.e. the Blue Plume. In Pap-I and (above all) in Sect. 2 
and 3.1 above we have shown that the BP population is a 
distinctive feature of the CMa system. Fig. 14 shows that 
this feature emerges clearly also from the GSC2.2 data: the 
well populated BP observed in the CMa direction (upper left 
panel of Fig. 14) has no counterpart in any of the other con- 
sidered fields. In particular the excess of BP stars (falling in 
the selection box reported in Fig. 14) in the CMa field with 

9 The same kind of problems forced us to limit the comparisons 
with GSC2.2 data to smaller fields with respect to the analysis of 
2MASS data. 




Figure 14. CMDs of the GSC2.2 selected fields described in 
Tab. 2. The overplotted selection box encloses the Blue Plume 
feature. 

respect to the WARP field is at the ~ 16a level, i.e. very 
large and significant. This finding provides additional sup- 
port to the conclusions of the previous section and indicates 
that there is at least a significant difference in the stellar 
populations hosted by the CMa system and by the Galactic 
warp. 

4.3 Kinematics 

M04 derived proper motions of M giants in the CMa over- 
density and found that the tangential velocity of these stars 
is compatible with that of Galactic disc stars at similar dis- 
tance. Here we simply note that this observation does not 
provide any conclusive argument on the nature of CMa since 
the uncertainties in the measured proper motions leave room 
for differences in tangential v elocity as large as ~ ± 70 km/s. 
In a dedicated contribution (Mar tin et alj|2004q| we show 
that, on the other hand, CMa stars display a systemic radial 
velocity significantly different from that of disc stars at the 
same distance, further supporting the idea that CMa is an 
independent stellar system. 



5 CLUSTERS AT THE SAME DISTANCE AS 
CMA 

In Pap-I it has been noted that there is also an ap- 
parent overdensity of open clus ters in the surround- 
ings of the CMa galaxy (see a lso iFrinchabov et aD 120041 : 
iForbes. Strader fc Brodiel l2004h . In the region of the sky 
considered in the present analysis we found useful photom- 
etry only for two of these clusters, namely Arp-Madore 2 
(AM-2) and Tombaugh 2 (Tom 2). 

According to the most recent studies llLed |l997; 
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lOrtolani. Bica fc Barbuv Ill995l) AM-2 is an old (age ~ 5 
Gyr) and relatively metal poor ([Fe/H] ~ —0.5) open clus- 
ter. Its position in the sky [(l,b)=(248.2°; -5. 8°)1 a nd its 
distance modulus ((m — M)o = 14.74±0.16|LeeJ[1997) place 
it fully within the CMa system, according to the distance 
estimates derived in the present analysis. Note also that the 
metallicity and age of this cluster are fully compatible with 
what found for the dominant population of the galaxy (see 
Sect. 2.5). Unfortunately the fiel d of view of the widest field 
photo metry currently available llOrtolani. Bica fc Barbuvl 

1995) is too small (5.5' x 5.5') to allow a comparative study 
with the population surrounding the cluster (supposedly 
mainly composed by CMa stars). The cluster population 
dominates the field (in the range of magnitudes covered by 
the availa ble CMP) at least out to the edg e of the field ob- 
served bv lOrtolani. Bica fc Barbuv I lll995f) though the den- 
sity profile of the cluster a ppears to be essentially flat for 
r > 1.5' (see, also lLee1IT997i) . 

The main characteristics of Tom 2 are quite similar to 
that of AM-2, but in this case the wider field of view (10.4' x 
10.4') of the best available photometry jKubiak et aljH992T) 
permits a more fruitful preliminary analysis, in the present 
context. Tom 2 has the same metallicity as AM -2 (\Fe/H] = 
— 0.5 ±0.1, from high resolution spectroscopy bv lBrown et alJ 

1996) and a similar age (~ 4 Gyr, the quoted age differ- 
ence with AM-2 is probably consistent with zero within 
the uncertainties). It is located at (l,b) = (232.9°; —6.9°) and 
its distance modulus ((m — M)o = 14.0 ± 0.3, according 
to lKubiak et alJll992ll is consistent with that of the CMa 
galaxy, within the errors. Also in this case the stellar popula- 
tion is similar at any distance from the cluster center (within 
the available field) though the profile is flat for r > 160" (see 
iKubiak et a!1ll992ri . 

In Fig. 15 the CMD of Tom 2 is compared with the CMD 
of F-XMM and with a synthetic CMD from the R03 model 
of the Tom 2 field. The observed and synthetic CMDs of 
Tom 2 have been shifted to the same reddening as F-XMM, 
correcting for the reddening differences indicated in the plot. 
The similarity of the Tom 2 population with that of CMa 
is striking and suggests an intimate connection between the 
cl uster and the galaxy, c onfir ming the hypothe s es adv anced 
bv IMartin et all ]2004al) and iFrinchabov et alJ ]2004l) . 

All of the results and discussion presented above sug- 
gests that AM-2 and Tom 2 are physically associated with 
the CMa system. A case strictly analogous to that of 
Tom 2, i.e. a cluster having a population indistinguish- 
able from the dominant population of the parent galaxy , 
is that of Ter 7 in the Sgr dSph fsee lMarconi et alJll997l) . 
This cluster is also similar to AM-2 and Tom 2 in age and 
metal content. For Tom 2 the hypothesis may also be ad- 
vanced that it is not a genuine cluster but a mere over- 
density within the Canis Major galaxy, of the kind that 
have been observed in the Sgr dSph l|lbata et alJ Il997t 
iBellazzi ni. Ferraro fc Buonanncj | l999al) and in the UMi 
dSph fs"ee TB^h^zzmT^t^dTl27ffl2bl and references therein). 
The issue can be definitely settled only with the observation 
of a significantly wider field surrounding the cluster. 
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Figure 15. The CMD of the open cluster Tombaugh 2 (middle 
panel, from Kubiak et al., 1992) is compared with the CMDs 
of the XMM filed (left panel) and of the synthetic CMD for an 
equivalent field from R03. The CMD of Tom 2 and the synthetic 
CMD have been shifted to the same reddening as F-XMM by 
applying the reddening differences indicated in the panels. 

6 SUMMARY AND CONCLUSIONS 

We have detected, for the first time, the Main Sequence of 
the main stellar population of the newly discovered galaxy 
in Canis Major (Pap-I), from photometry of a wide field 
located at ~ 4.2° from its center (the XMM field). From the 
analysis of the CMD of the XMM field we conclude that: 

• The stellar population observed in the XMM field has 
no counterpart in the R03 Galactic model. 

• The distance modulus of CMa is (m-M) ~ 14.5 ±0.3, 
where the reported uncertainty mainly reflects the uncer- 
tainty in the actual average metallicity of the system, that 
we assumed to lie in the range — 0.3 ^ ]MJH] ^ —0.7, 
according to iMartin et alJ j2004al) : ICrand {2003). The cor- 
responding distance is Dq ~ 8.0 ± 1.2 kpc. 

• The age of the main population of CMa lies in the range 
4 < age< 10 Gyr, slightly depending on the assumed metal- 
licity. 

• There is a sparse population of younger stars (age< 1 — 2 
Gyr) that is also associated with the CMa system (Blue 
Plume) . 

Tentative detections of stars related to the CMa/Ring 
system are presented for the fields surrounding the open 
clusters NGC 2477, Tombaugh 1 and Berkeley 33. The cases 
of two open clusters (AM 2 and Tombaugh 2) that may be 
physically associated with the CMa system have also been 
discussed. The results presented in this study provide in- 
dependen t confirmation full s upport and complement the 
results of IMartin et al] <l2004al) . that described for the first 
time the CMa ste llar system. 

The claim bv lMomanv et al] J2004I) that the CMa over- 
density could be entirely due to the South/North asymme- 
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try produced by the Galactic warp has been demonstrated 
to lack the essential observational basis. Using the F-XMM 
data we have shown that the warp model encapsulated in 
the R03 model (used by M04 to substantiate their claim) 
fails to reproduce both the star counts and CMD morphol- 
ogy observed in this field, located at ~ 4.2° from the center 
of the CMa system. Using Red Clump and M giant stars 
from 2MASS we have shown that (a) the CMa overden- 
sity is much larger than that produced by the warp, and 
(b) it is found only at a well defined distance not compati- 
ble with the warp, fully confirming the results of Pap-I and 
demonstrating that if all the available information is taken 
into account the CMa overdensity clearly emerges above the 
"noise" provided by the Galactic warp. Using GSC2.2 data 
we demonstrated that the Blue Plume population is a dis- 
tinctive characteristic of the CMa system that has no com- 
parable counterpart neither in Northern Control Fields or 
in fields sampling the real warp. 

A better characterization of the stellar content and star 
formation history of the CMa system must wait for dedi- 
cated photometry and spectroscopy fully suited for the study 
of this difficult target, immersed in very contaminated and 
extincted regions of the sky. Here we tried to obtain the 
best possible characterization of the system from existing 
and publicly available data, as an immediate follow-up of 
the discovery. 

Future searches for CMa stars in the background of 
other targets may prove extremely useful in characteriz- 
ing the spatial distribution of the galaxy and, in particu- 
lar, of the associated Ring. The relevant change of view is 
that what may be perceived as an undesired contaminant of 
a given observation may hide more astrophysical informa- 
tion than the target itself. This is true either if the back- 
ground object is an unknown stellar system (as CMa or 
Sgr) either if it is a poorly characterized Ga lactic compo- 
nent, as the stellar warp of the Milky Way feinnevlll992t 
iKuiiken fc Garcia- Ruid lioflF). The lessons provided by the 
Sgr dSph (and Stream) and by the CMa galaxy (and Ring) 
suggests to all of us a renewed attitude in the analysis of 
astronomical data, to avoid missing relevant pieces of infor- 
mation hidden in the background of our main targets. 
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